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ABSTRACT 

We present a comprehensive statistical analysis of Swift X-ray light-curves of Gamma- 
Ray Bursts (GRBs) collecting data from more than 650 GRBs discovered by Swift and 
other facilities. The unprecedented sample size allows us to constrain the rest-frame 
X-ray properties of GRBs from a statistical perspective, with particular reference to 
intrinsic time scales and the energetics of the different light-curve phases in a common 
rest-frame 0.3-30 keV energy band. Temporal variability episodes are also studied and 
their properties constrained. Two fundamental questions drive this effort: i) Does the 
Jv> ' X-ray emission retain any kind of "memory" of the prompt 7-ray phase? ii) Where 

'j_j I is the dividing line between long and short GRB X-ray properties? We show that 

C^ . short GRBs decay faster, are less luminous and less energetic than long GRBs in 

the X-rays, but are interestingly characterized by similar intrinsic absorption. We 
furthermore reveal the existence of a number of statistically significant relations that 
link the X-ray to prompt 7-ray parameters in long GRBs; short GRBs are outliers 
of the majority of these 2-parameter relations. However and more importantly, we 
report on the existence of a universal 3-parameter scaling that links the X-ray and 

the 7-ray energy to the prompt spectral peak energy of both long and short GRBs: 
T? „ 77i.oofo.O6 /pO.eoio.io 

iiX.iso OC ii^ i^^ /iSp]^ 

Key words: gamma-ray: bursts - radiation mechanism: non-thermal -X-rays 



1 INTRODUCTION (GRBs): the X-ray afterglow . The sta ndard afterglow the- 

.,. „ ,. . „ , rr—. — ■ n u^^ J n , ory (Mes zaros fc Re cs 1997; Sari et al. 1998) predicts that 

In 7 years or operation, bwift HGehrels et al.ll2004l l has rev- .^ ' . . . . ,, . , , • r 1 x- • i- 

■^ r- J J 1 a/ A-ray emission arises rrom the interaction 01 a relativistic 

olutionized our understandina; of the X-ray emission that . „ -^1^1 1 • ^ i- 1 i- ^ ^i r 

"■^ outnow with the ambient medium, leading to the rorma- 
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tion of a blast wave. In this context, short and long GRBs 
would naturally show similar afterglows (since the emis- 
sion would be sensitive to the energy budget of the rela- 
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tivistic outflow but otherwise keep no memory of the ori- 
gin of the outflow) if the properties of the environment of 
the two classes are also similar. This is difficult to recon- 



cile with the massive star (long GRBs, see e. g. IWooslev 
1993 ) vs. compact binary (short GRBs see e.g. IPaczvnski 



1986f ) progenitor systems, which would instead suggest a 



wind density profile (i.e. oc r~ ) around long GRBs and 
an ISM (i.e. oc r°) ambient density for short bursts. Con- 
trary to expectations, observations are often consistent with 
a constant densi ty environment (ISM ) also in the case of 
long GRBs (e.g. iRacusin et al.l 120091) . Moreover, the stan- 
dard afterglow theory fails to explain the presence of long 
(~ 10'' s) phases characterized by very mild decays (the 
so-called plateaus or shallow decay phases) i n the X-ray 
light- curve of many GRBs (e.g. GRB 060729, iGrupe et all 
[20071); it cannot account for abrupt dro ps of e mission ob- 
served in some G RBs (e.g. GRB 070110. iTroia et al. 2007, 
iLvons et al.l |2010|'1 and has serious d i fficulties explaining 
the X-ray flares jChiricarini et al.ll2007l . iFalcone et al.ll2007l : 
IChincarini et al.ll2010l and references therein). 

As a result, a number of alternative models have been 
proposed. They basically divide into two classes: accretion 
onto a newly born black hole which directly powers the ob- 
served X-ray light-curve (JKumar et al.ll2008|) or power f rom 
a rapidly rotating magnetar (see e.g. lMetzger et al.ll201ll and 
references therein). In sharp contrast to the standard after- 
glow theory, those models directly relate the properties of 
the observed X-ray light-curves to the GRB central engine. 

With this work we improve our understanding of the 
X-ray emission of long and short GRBs through a homo- 
geneous analysis of a sample of more than 650 GRBs ob- 
served by Swift-XKT (|Burrows et al.l 120051 ). We ask: what 
is the typical amount of energy released during the different 
X-ray light-curve phases? Does the X-ray emission retain 
any kind of "memory" of the prompt 7-ray phase? GRBs 
are traditionally classified into long and short according to 
their prompt 7-ray properties: do short GRBs show a dis- 
tinct behaviour in the X-rays as well? Is it possible to find 
a universal (i.e. common to long and short GRBs) scaling 
that involves prompt and X-ray properties? This set of still 
open questions constitutes the major reason to undertake 
the present investigation. 

Previous attempts mainly concentrated on observer 
frame properties and tried to understand to what extent 
the observations could b e reconciled with the standard for- 
ward shock model (e.g jO'Brien 6^111120061: [Butlej |2007|; 



Butler fc Kocevskill2007l: IWillingalc et al 



2007 



20071 . I2OO8I : lEvans et al] 12009: Racu sin et al 



. Liang et al.l 
20091 . l201lh . 



This effort led to the identification of serious difficulties of 
the standard picture. We build on previous results and adopt 
here a different approach: instead of comparing observations 
with a particular physical model, we take advantage of the 
large sample size and look for correlations between the X- 
ray and 7-ray properties of GRBs any physical model will 
have to explain. We complement previous studies with: 

i. Homogeneous analysis of GRBs in a common rest frame 
energy band (0.3-30 keV). 

ii. Statistics and properties of the temporal variability su- 
perimposed on the smooth X-ray decay. 

iii. Comparative study of long vs. short GRB X-ray af- 
terglows. 



iv. Study of the prompt 7-ray vs. X-ray connection: no- 
tably, we report on the existence of a universal scaling in- 
volving prompt and X-ray parameters. 

Hereafter, we will refer to the X-ray signal recorded by 
Swift-XKT after a GRB trigger as "X-ray light-curve" (LC) 
and explicitly do not use the word "afterglow" to avoid con- 
fusion ("afterglow" refers to the standard interpretation). 
Uncertainties are given at 68% confidence level (c.l.) unless 
explicitly mentioned. Standard cosmological quantities have 
been adopted: -ffo ~ 70kms~^ Mpc~', Qa = 0.7, SIm = 0.3. 
The results from our analysis are publicly availablelj 



2 DATA ANALYSIS 

We select GRBs observed by Swift-XKT from the beginning 
of science operations in December 2004 through the end of 
2010. The starting sample includes 658 GRBs; 36 belong t o 
the class of short GRBs. Following [Margutti etljl (|2011bh . 
the short or long nature of each event is established using 
the combined information from the duration, hardness and 
spectral lag of its prompt 7-ray emission: a prompt 7-ray 
duration Tgo < 2 s coupled to a hard 7-ray emission with 
photon index F < 1.5 and a negligible spectral lag are con- 
sidered indicative of a short GRB nature. 

For each GRB, the data reduction comprises four steps: 

• Extraction of count-rate light-curves (LC hereafter) in 
the 0.3-10 keV energy band. 

• Time-resolved spectral analysis. 

• Flux calibration in the observer frame 0.3-10 keV en- 
ergy band and luminosity calibration in the rest-frame 0.3-30 
keV energy range for the sub-sample of GRBs with known 
redshift. 

• LC fitting. 

The Swift-XRT data have been analyzed using the lat- 
est version of the HEASOFT package available at the time 
of the analysis (v. 6.10). For each GRB we started from 
calibrated event lists and sky images as distributed by the 
HEASARC archivqj. The following analysis made extensive 
use of the XRTDAS software package. The additional au- 
tomated processing was performed via custom IDL scri pts. 
Detail s on the procedure followed can be found in iMarguttil 
(|2009l )[^ Here we note that: 

(i) While agreeing on the major steps, the data extraction 
a dopted here is slight ly diff erent from the meth ods presented 
in [Evans et al.l (|2007l ) and [Evans et al.1 l[2009l ): however, af- 
ter comparing every single LC obtained with the two tech- 
niques, we find that the methods lead to consistent results. 

(ii) Our flux and luminosity LC calibration is based on 
a time-resolved spectral analysis which is able to capture 
the spectral evolution of the source with time. Uncertainties 
arising from the spectral analysis have also been propagated 
into the final flux and luminosity LCs (this is essential to 
compute the significance of positive temporal fiuctuations 
superimposed on the smoothly decaying LC, see Sec. l2.ip . 



^ A demo version of the website is currently available at 
'http://www.grbtac.org/xrt_denio/GRB060312Afterglow.h tnil] 
^ http://heasarc.gsfc.nasa.gov/cgi-bin/W3Browse/swift.pl 
3 Retrievable from h ttp://hdl.handle.net/10281/7465l 
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Figure 1. Cartoon showing tiie different LC models used in tiiis 
work. Botti axes use logaritiimic units. 



(ill) The intrinsic neutral hydrogen absorbing column NHhg 
of each GRB was computed extracting a spectrum in the 
widest interval of time with no apparent spectral evolution. 
The best-fitting NHhg was used as frozen parameter in the 
time-resolved spectral analysis above. The Galactic contri- 
bution was froz en to the value i n the direction of the burst 
as computed bv lKalberla et al.l (|2005l ). 

Since we corrected for the Galactic and intrinsic ab- 
sorption, the final results are unabsorbed 0.3-10 keV (ob- 
server frame) flux LCs. For the sub-sample of GRBs with 
known redshift we furthermore extracted unabsorbed lumi- 
nosity LCs in the 0.3-30 keV (rest-frame) energy band (ex- 
trapolating the best-fitting power-law spectrum). We con- 
servatively use only z derived from optical spectroscopy and 
photometric redshifts for which potential sources of degen- 
eracy (e.g. dust extinction) can be ruled out with high con- 
fidence. The complete list of redshifts used is reported in 
Table IBTI (175 GRBs in our sample have redshift). 



2.1 Light-curve fitting 

The X-ray LCs of GRBs for t > 60s consist of smoothly 
decaying power-laws or broken power-laws with X-ray 
flares superimposed. Here we concentrate on the underly- 
ing smooth component. 

We considered only GRBs whose statistics were good 




Obs Time since BAT trigger (s) 



Figure 2. 0.3-10 keV unabsorbed flux LC of GRB 060312 witfi 
best-fitting model superimposed (blue solid line). The red crosses 
mark the data points identified as "excesses" during the semi- 
automatic fitting routine and subsequently removed from the fit 
to obtain the best-fitting power-law plus broken power-law model. 
Vertical dashed lines: best estimates of the break times obtained 
as described in the main text. 



Table 1. Number of GRBs per LC type. For each type, the num- 
ber of GRBs with complete LCs (C-GRBs, defined as promptly 
repointed GRBs Uep < 300 s whose fading was followed up down 
to a factor ~ 5 — 10 from the background limit) and with detected 
positive fiuctuations (F-GRBs) with respect to the best-fitting 
smooth decay is also reported. GRBs are classified as either C or 
U-like, and either F or N-like. U-like bursts have truncated LCs, 
while N-like GRBs show no evidence for fiares. This classification 
refers to the 0.3-10 keV (observer frame) LCs. 



Light-curve types 







la 


lb 


Ila 


lib 


III 


Total Number of GRBs 


114 


89 


61 


133 


18 


22 


C-GRBs 


42 


61 


53 


121 


17 


22 


F-GRBs 


23 


16 


24 


48 


8 


10 



enough to allow us to extract a spectrum to convert their 
count-rate LCs into flux LCs (total of 437 GRBs out of 658) . 
We first fitted the entire sample of flux LCs in the 0.3-10 keV 
(observer frame) energy band. We then focussed on the sub- 
sample of GRBs with redshift and performed a second fit 
using the LCs in the common 0.3-30 keV (rest-frame). 

Our semi-automatic fitting routine is based on the 
x'^ statistic and closel y follows the procedure outlined in 
iMargutti et al.l (12011^ ). We fit the foUowing models. Defin- 
ing: 



/(iVi,ai,i) = 7Vit- 



(1) 



(where A^i and ai are the normalization and the slope of 
the power-law, respectively) and 



(2) 



g{N2,a2,a3,ti,s,t) = N2[ ( — 



+ U 



(where N2 is the normalization, 02 and as are the slopes 
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of the broken power-law; i(, is the break time while s is the 
smoothing parameter), the fitting models can be written as: 



Short GRBs 



• Simple power-law (model 0): 

F = /(7Vi,ai,i) 



(3) 



• Smoothed broken power-law (model la and lb for s < 
and s > 0, respectively): 



F = g(Ni,ai,a2,tti,si,t) 



(4) 



• Smoothed broken power-law plus initial (model Ila) 
power-law decay: 



F = f{Ni,ai,t) + g{N2,a2,a3,tb2,si,t) 
or final (model lib) power-law decay: 

F = g(Ni,ai,a2,tti,si,t) + f(N2,a-i,t) 



(5) 
(6) 



• Double smoothly joined broken power-laws (model III): 

F = g{Ni,ai,a2,tt-i_,Si,t) + g{N2, 03,04,^3, S2,t) (7) 

The model number follows from the number of break times. 
For model Ila (lib, III) the first (second) break time is 
defined as the time when the second (first) component out- 
shines the first (second) component. Figure [T] illustrates 
the different models, while Fig. [2] shows the result for 
GRB 060312 taken as an example: in this case the semi- 
automatic procedure identified two episodes of emission in 
excess of the smooth decay (model Ila). The number of 
GRBs for each LC type are listed in Table [l] We refer to 
the LCs as "type" 0, la, lb, Ila, lib and type III GRBs in 
the following. 

The best-fitting parameters together with their uncer- 
tainties and associated covariance matrix were then used to 
derive the 0.3-10 keV (observer-frame) fiuence of the entire 
LC, from the S'uii/i-XRT re-pointing time to the end of the 
observation. No temporal extrapolation was applied at this 
stage. Note that the contribution from significant positive 
fluctuations has not been included. The fluence of the dif- 
ferent LC phases as defined by the temporal breaks was also 
calculated (Fig. [l|. Results are listed in Table iBll and Table 
IB3I We then followed the very same procedure to fit the 0.3- 
30 keV (rest-frame) LCs: Table lB5l reports the energetics in 
this energy range. 

The list of LC points flagged as "excesses" during the 
fltting procedure (e.g. red crosses of Fig. [2| constituted for 
each GRB the starting sample to look for significant posi- 
tive fluctuations with respect to the best fit. The information 
contained in the covariance matrix was used to derive the un- 
certainties associated with the residuals with respect to the 
best-fit (residuals were at this stage calculated on the entire 
LC). We first selected positive fluctuations with a minimum 
la significance. We furthermore require the positive fluctua- 
tions to show a rise plus decay structure: this procedure au- 
tomatically excluded single data-points scattering from the 
best fit. GRBs showing (not showing) such structures were 
fiagged as "F" ("N") in Table [BT]|B5l GRB 060312 in Fig. 
O with two rising and decaying structures superimposed on 
the smooth decay qualifies as "F"-event. The fiuence (energy 
for known z) of those excesses was calculated by simply inte- 
grating the flux of each LC bin over the bin duration (after 
subtracting the contribution from the underlying smoothly 
decaying emission). Errors were propagated accordingly and 



050724 051221A 051227 060313 061006 
061201 070714B 070724A 070809 071227 
080123 080503 080919 090510 090515 
090607 100117A 100816A 101219A 

Table 2. List of 19 short GRBs with complete LCs. GRBs 
with detected tem porally extended emission are in boldface 
JNorris et aIJl2011^ . 



can be used to quantify how significant is the presence of 
emission in addition to the smooth power-law decay in each 
GRB. The fluence (energy) of positive fiuctuations detected 
during the different LC phases (e.g. steep-decay, plateau, 
normal-decay, etc.) was also derived and listed in Table lB3l 
(0.3-10 keV, observer frame) and Table |B5] (0.3-30 keV, rest 
frame). For simplicity, in the following we will use the word 
"flare" to refer to statistically signiflcant positive fluctua- 
tions detected on top of the smoothly decaying component, 
being however aware that different kinds of variability pos- 
sibly contribute to the detected "flaring activity" . 



3 LONG VS. SHORT GRBS PROPERTIES 

The analysis above reduces the X-ray LC of GRBs to a set 
of measured parameters: temporal slopes; break times; to- 
tal isotropic energy (fluence) and energy (fiuence) associated 
with the different LC phases; flare energy (fluence); spectral 
photon index temporal evolution; intrinsic neutral hydrogen 
absorption. This constitutes an unprecedented set of infor- 
mation homogeneously obtained on the largest sample of 
GRBs to date and represents the natural sample to look for 
correlations among the parameters. 

While we report the best-fitting parameters of the en- 
tire sample (Tables IB1IIB5|I , in the following we restrict our 
analysis to GRBs with "complete" LCs, deflned as those 
GRBs re-pointed by XRT at frep < 300 s and for which we 
were able to follow the fading of the XRT flux down to a 
factor ~ 5 — 10 from the background limit (or, equivalently, 
tend > 4 X lO'^ s). These GRBs are flagged as "C" in Tables 
[BT]|B5l The number of "C"-hke GRBs per LC morphologi- 
cal type is reported in Table [T] "U"-like GRBs have instead 
truncated LCs. Short GRBs with complete LCs are listed in 
Tab le El Short GRBs w ith extended emission are in boldface 
(see iNorris et allboill l . 



3.1 Median X-ray light-curve of long and short 
GRBs 

We select the sub-sample of C-like LCs of GRBs with 
redshift observed in the common rest-frame time interval 
10^ - lO'^ s and 10^ - 10* s for long and short GRBs, respec- 
tively. These criteria resulted in a sample of 79 long GRBs 
and 9 short GRBs (Fig. |3}. We consider here 0.3-30 keV 
(rest-frame) luminosity LCs. 

We combined the best-fltting profiles of Sec. l2.1l to pro- 
duce a median luminosity LC of a long GRB. The result is 
shown in Fig. the median luminosity LC roughly decays 
as ~ t^^ (with a milder decay ~ t~°'^ for (0.5 < t < 4) ks 
and a steeper ~ t~^'^ decay after ^ 4ks). With the possible 
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Figure 3. Best fitting profiles of 9 C-like siiort GRBs witii rest- 
frame time coverage 10^ — lO^s (colour lines) superimposed on the 
sample of long GRBs (grey lines). Black (blue) solid line: median 
Log{L) light-curve for short (long) GRBs; dashed lines mark the 
la dispersion. The LCs have been calibrated in the rest-frame 
0.3-30 keV energy band. 



Figure 4. 0.3-30 keV (rest-frame) X-ray energy as a function of 
redshift. Black (red) points; long (short) GRBs. Blue dashed line: 
empirically derived detectability threshold. For z > 2 we are not 
sensitive to GRBs with Ex iao S lO^'^ erg. For z > 1 there is 
no evidence for an evolution of the upper bound of Ex iso with 
redshift. 



exception of the shallower section, this is in rough agree- 
ment with the prediction of the standar d afterglow theorjO 
IJMeszaros fc Reejl 19971 : ISari et al.l[l99i ) . 

The decay of the median LC is steeper for short (oc 
t~^'^ ) than for long GRBs (oc t~^) in the rest-frame time 
interval 10^ — 10* s; short GRB X-ray LCs are on average 
less luminous by a factor ~ 10 — 30 than long GRBs X-ray 
LCs. This conclusion holds also considering long GRBs in 
the same redshift bin. However, Fig.[3] clearly shows that th e 
two samples slightly overlap (see also ICehrels et al.l 120081 ) . 
The steeper decay that characterizes short GRBs causes a 
progressive shift of their luminosity distribution towards the 
low end of the long GRB distribution. 



3.2 Energetics of long and short GRBs 

Table [3] reports the analysis of the parameter distributions 
derived from the LC fitting of Sect. [2l A complete list of 
symbols can be found in Appendix [X] The observed Ex,iso 
distribution peaks at ~ 7 x 10^^ erg, typically represent- 
ing ~ 7% of the 1 - 10* keV (rest-frame) E^^i^o- Figure g] 
shows that we are not sensitive to the population of bursts 
with -Ex, iso < lO'""^ erg for z > 2 (so that the low energy 
tail of the Ex.iso distribution is currently under-sampled). 
This is likely a non-detectability zone, consequence of the 
Ex, iso oc i5!J;fso of Sec. 21 For z > 1 there is no evidence for an 
evolution of the upper bound of Ex, iso with redshift, which 
may suggest that ~ 10^'^ erg is a physical boundary to the 
Ex, iso distribution (the record holder is GRB 080721 with 
Ex, iso ^ 10^^ erg)- lu this respect we note that maximum 
budget Emax ~ 10'''^ ergj is predicted by magnetar mod- 



4 A steepening to ~ 1.5 — 2 is predic ted after the jet-break tim e 
if the outflow is coUimated into a jet l lRhoadall999l : ISarilll999r) . 
^ It is not a given that GRB 080721 violates this limit: -Ex, iso 
represents the isotropic equivalent X-ray energy, an overestimate 
to the true value if the emission -as we believe- is beamed. 



els (|Usovlll992l ). The same pattern is followed by the flare 
energy E^^: for 2 > 2 we are not sensitive to E^^ < 10^'' 
erg. 

Observations suggest tha t the GRB X-ray LCs consist 
of two distinct phases (see e.g. lWillingale et al.ll2007l ): a first 
steep dec ay phase tightly connected to the prom pt 7-ray 
emission (|Tagliaferri et al.ll2005l : iGoad et al.l 1200^ ): and a 
second phase characterized by a flattening of t he LC (with 
limit ed evidence for spectral evolution, see e.g. iLiang et al.l 
120071 ) followed by a "normal decay" phase. Type Ila GRBs 
(Fig. [1} clearly show the presence of both components, with 
energy Ei,x = Ei and E2,x ~ E2 + E3, respectively; for type 
la LCs, the lack of spectral evolution and the typically mild 
slope ai resembling 02 of type Ila lead us to identify E2,x ~ 
E1 + E2; type lb GRBs show strong spectral evolution during 
the flrst LC segment: this together with the transition to a 
milder decay at i(„i lead us to define Ei,x = Ei and E2,x ~ 
E2; in type lib LCs the spectral and temporal properties 
of the first segment (with slope ai) strongly suggest that 
Swift-yiRT caught the end of the prompt emission in the 
X-rays: we therefore define Ei,x = Ei + E2 and E2,x = 
E3; the same is true for type III GRBs: in this case we 
define Ei,x — Ei + E2, E2,x — Eg + E4. The two phases 
release comparable energy (see Table[3]), with Ei,x and E2,x 
peaking at ~ 1.1 x 10^^ erg and ~ 4 x 10^^ erg, respectively. 

In each distribution, short GRBs populate the low en- 
ergy tail: Ex'^i°o' ~ 10^" erg, which is approximately 2 orders 
of magnitude less than a typical long GRB. Figure |4] also 
shows that short GRBs are less energetic than long GRBs 
in the same redshift bin. A systematic difference between 
the -still poorly constrained- jet opening angles of long and 
short GRBs, with short GRBs being less coUimated than 
long GRBs, could in principle mitigate this energy gap. If 
we compare the energy released during the two phases sep- 
arately (i.e. early steep decline vs. plateau plus subsequent 
decay), we find an indication that short GRBs are more en- 
ergetically deficient during the second phase then in the first 
phase, i.e. Elt^VE 



long 
X.2 



0.014 and E^^V E[ 



long 



0.054. 
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Table 3. Characteristic quantities describing the parameter distributions (number of elements (#), mean (rra), median (Af), standard 
deviation (SD), skewness (SK)) and best-fitting values from a Gaussian fit (mean (fj,), standard deviation (tr), normalization (N)). 



Fluences (S) are given in 10 erg cm , energies (E) in 10 erg, fluxes (F) in 10 



s cm , luminosities (L) in 10 "ergs 



times (t) in s, hydrogen column densities (NH) in lO^'^ cm~-^. Note that logarithmic (linear) units have been used in the upper (lower) 
half of the Table. We refer the reader to Appendix |X] for the exact definition of the parameters listed below. X-ray energies, luminosities 
and intrinsic times have been computed in the rest-frame 0.3-30 keV energy band. All the other X-ray quantities refer to the 0.3-10 
(observer frame) band. 





# 


m 


M 


SD 


SK 


/^ 


(T 


N 


Log{S^) 


386 


0.17 


0.18 


0.61 


-0.07 


0.17 ±0.05 


0.60 ±0.04 


105.1 ±7.7 


Log{Tgo) 


334 


1.58 


1.73 


0.65 


-1.18 


1.67 ±0.07 


0.62 ±0.06 


93.9 ±9.4 


Lo9(El,^-^^^) 


151 


2.06 


2.37 


0.93 


-1.12 


2.26 ±0.14 


1.02 ±0.1 


91.6 ±9. 


Log{E^^iso) 


78 


2.88 


3.01 


0.91 


-0.89 


3.01 ±0.14 


0.85 ±0.12 


38.9 ±5.7 


Log{Epi^) 


78 


2.64 


2.71 


0.52 


-0.76 


2.56 ±0.1 


0.47 ±0.069 


23.6 ±4.7 


iog(Lpk,iso) 


85 


2.43 


2.51 


1.00 


-1.90 


2.51 ±0.19 


0.91 ±0.12 


61.0 ±9.3 


iog(TgRP) 


138 


1.18 


1.33 


0.59 


-0.65 


1.26 ±0.11 


0.67 ±0.09 


53.8 ±7.4 


Lo9(NHhg) 


161 


21.6 


21.8 


1.21 


-3.54 


21.9 ±0.1 


0.62 ±0.09 


42.5 ±7.3 


Log{Sx) 


316 


-0.38 


-0.42 


0.62 


0.23 


-0.46 ±0.05 


0.57 ±0.04 


68.6 ±5.3 


Log(Sl^) 


115 


-0.81 


-0.76 


0.82 


-0.38 


-0.76 ±0.12 


0.89 ±0.12 


51.3 ±6.0 


LogiS^^y,) 


211 


-0.90 


-0.92 


0.80 


0.12 


-0.89 ±0.10 


0.88 ±0.06 


97.2 ±7.9 


Log{S2,x) 


316 


-0.58 


-0.61 


0.63 


0.14 


-0.57 ±0.07 


0.70 ±0.06 


102.2 ±8.5 


Log{Sf\) 


62 


-0.75 


-0.69 


0.73 


-0.38 


- 


- 


- 


Log(Sl^) 


71 


-1.08 


-0.98 


0.90 


0.03 


-1.08 ±0.17 


0.95 ±0.16 


36.7 ±5.6 


Log{Ex,iso) 


126 


1.67 


1.84 


0.81 


-0.67 


1.82 ±0.08 


0.88 ±0.08 


31.1 ±2.5 


Log{El^) 


59 


1.25 


1.40 


0.97 


-0.64 


- 


- 


- 


LogiEi^x) 


86 


1.00 


1.04 


0.92 


-0.28 


1.10 ±0.11 


0.94 ±0.08 


40.7 ±3.6 


Log(E2,x) 


126 


1.45 


1.63 


0.92 


-0.94 


1.63 ±0.11 


0.82 ±0.10 


63.2 ±7.0 


Log(Ef\) 


35 


1.13 


1.38 


1.00 


-0.78 


- 


- 


- 


Log{E^^) 


38 


1.04 


1.14 


0.98 


0.01 


- 


- 


- 


Log(ti) 


155 


2.66 


2.56 


0.48 


1.06 


- 


- 


- 


Log{t{) 


155 


3.94 


3.93 


0.73 


0.19 


3.93 ±0.14 


0.8 ±0.12 


59.2 ±8.9 


Log(FO 


155 


-4.23 


-4.26 


0.83 


-0.14 


-4.18 ±0.11 


0.89 ±0.12 


61.1 ±6.6 


Log{Fi) 


155 


-5.01 


-4.94 


0.80 


-0.12 


-4.99 ±0.11 


0.80 ±0.08 


63.4 ±6.9 


Logitf^) 


62 


2.13 


2.03 


0.62 


1.27 


- 


- 


- 


Logitf^) 


62 


3.58 


3.48 


0.74 


0.29 


3.53 ±0.11 


0.75 ±0.09 


23.0 ±2.9 


Log{Li) 


62 


0.54 


0.73 


1.25 


-1.07 


0.79 ±0.19 


1.26 ±0.16 


46.6 ±5.5 


Log(L{) 


62 


-0.47 


-0.19 


1.19 


-1.19 


0.04 ±0.33 


1.41 ±0.29 


38.2 ±6.8 


cist 


213 


3.96 


3.56 


2.34 


3.95 


3.22 ±0.51 


2.34 ±0.35 


238.0 ±35.0 


"sh 


155 


-0.16 


0.18 


1.23 


-4.06 


0.27 ±0.14 


0.52 ±0.12 


61.0 ±13.0 


an 


204 


1.59 


1.38 


1.04 


8.30 


1.34 ±0.13 


0.49 ±0.11 


93.0 ±19.0 



This argues against a beaming related explanation, since the 
jet opening angles of long and short GRBs are expected to 
be more similar at late than at early times. Short GRB light- 
curves decay faster than long GRBs in the X-rays, typically 
resulting in shorter observations {t^nd ^ 10"* s vs. t^„d ^ 10^ 
rest-frame): however, using the average L oc t~^'^ scaling 
above, we find E^^^^'^t < 10* s) ~ E^xTi^ < 10^ s). Thus 
the relatively lower measured energy of the later LC phase 
in short GRBs compared to long GRBs is not due to the 
shorter observations. 



3.3 Intrinsic neutral hydrogen absorption in long 
and short GRBs 

The distribution of the intrinsic neutral hydrogen columns 
NHhg is portrayed in Fig. [S] panel (c). The distribution 
of measured NHhg is found to have an average valuqj 



^ Solar abundances are used to determine the best- fitting NHhg ■ 



m = 10 cm~ . When modeled with a lognormal distri- 
bution, the best fitting mean and standard deviation are: 
fj, — 21.9 ± 0.1, a = 0.6 ± 0.1, i n agre ement with the esti- 
mates by ICampana et al.l (|2010| . |2012[) obtained on smalle r 
samples. However, differently from ICampana et al.l (|2010|): 
(i) our sample contains a larger number of GRBs for which 
no evidence for intrinsic absorption was found (upper limits 
in Fig. O; (ii) we find evidence for a larger population of 
highly absorbed (NHhg > lO^^cm"^) GRBs at low redshift 
(z<2). 

A trend for increasing NHhg with redshift is apparent in 
Fig- El panel (a): however, our sensitivity to small amounts 
of intrinsic absorption decreases with increasing redshift due 
to the fixed XRT band-pass, which explains the higher per- 
centage of upper limits in the 4 to 6 redshift interval, and 
is at least partially responsible for the observed trend. The 
sample is furthermore redshift selected, which implies a bias 
against highly extinguished GRBs. The severity of this bias 
is possibly redshift dependent, with a dependence which is 
difficult to quantify. 

Short GRBs map the low end of the NHhg distribution. 
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Figure 5. Left panel: (a): Intrinsic neutral hydrogen absorption 
vs. redshift for long and short GRBs (grey and red dots, respec- 
tively). 90% upper limits are marked with arrows. Median NHjjG 
values in different redshift bins are indicated with filled trian- 
gles: for each bin, the error bars span the Icr NHhg dispersion. 
(b): redshift distribution of the sample of long (yellow) and short 
(red) GRBs. (c): Intrinsic neutral hydrogen distribution for long 
(black line) and short (red line) GRBs. The dashed histogram 
includes upper limits. In panels (b) and (c) a dashed line indi- 
cates the median value for the entire distributions {{z) = 1.82, 



(NHhg> = 102 



')■ 



with an average absorption NHIJq'^* = lO^^'* cni"^ (mean 
of the logarithm of NHIJq'^*). Their properties are however 
consistent with the intrinsic absorption of long GRBs in the 
same redshift bin. A KS-test comparing the NHhg distribu- 
tion of long and short GRBs with < 2: < 1 reveals that 
there is no evidence for long GRBs to show higher NHhg 
when compared to short GRBs in the same redshift bin (KS 
probability of 34% ). A possibility is that we missed the pop- 
ulation of long GRBs with even higher NHhg at low z (see 
above). This would imply that GRBs with low optical ex- 
tinction but high NHhg are typical of the high-redshift uni- 
verse, only (Wats on fc Jakob sson 201^). We conclude that 
using the available data, caution must be therefore used to 
interpret the long GRB NHhg distribution as aproof of their 
association to star formation (|Campana et al.ll2010l ) unless 
this association is meant to be extended to short GRBs as 
well. 



4 PARAMETER CORRELATIONS 

Here we proceed to look for 2-parameter correlations in- 
volving both X-ray and 7-ray properties. From a blind anal- 
ysis we found 199 statistically significant correlations (out 
of 946). We focus on the physically interesting, correlations. 
The significance of each correlation is estimated using the 
R-index r, the Spearman rank p and Kendall coefficient K 
(Table|3)). Only correlations for which the chance probability 
associated with at least one of the test statistics is < 10~^ 
have been listed. As a general note: 



• No significant correlation is found to involve the rest 
frame Tgo, the intrinsic NHhg or the LC temporal slopes; 

• We re-scaled the LC temporal breaks fb by the Tgo, 
adding new parameters to our list: y = t^/Tgo- However, 
the use of re-scaled properties did not improve any of our 
correlations and are therefore not included in the following 
discussion. 

The correlation coefficients and the best-fitting power- 
law parameters of each correlation are listed in Table U 
Our best-fitting p rocedure accounts for the sample variance 
(JD'Agostinill2005l ). 



4-0.1 The link between the X-ray and prompt y-ray energy 

Fig- El (left panel) shows that Ex,iso is directly linked to the 
isotropic energy released in 7-rays during t he prompt emis- 
sion i ?.T,i3o. A similar result was found by IWillingale et al.l 
(I 2OO7I ) on a smaller sample of GRBs. Here we show for 
the first time how short GRBs compare to long GRBs: no- 
tably, all short GRBs but GRB 050724 are outliers of the 
long GRB relation, with -Ex.iso for short GRBs a factor 
^ 50 below that for long GRBs and having large disper- 
sion (the i?x,iso distributions are almost distinct for long 
and short GRB s, as shown in Fig, lit . A clear exception 
is GRB 050724 (JBarthelmv et al.l bOoH : iGrupe et al.l I2OO6I ) 
which had a bright and long-lived X-ray afterglow with a 



powerful late time re-brightening (iBernardini et al.l I2OIII : 
ICampana et al.ll2006l : iMalesani et al.ll2007l ). This difference 
may be understood in terms of a different radiative efficiency 
rj-y (where rj-y = E~,/{E^-\-Ek), being Ek the outfiow kinetic 
energy) during the prompt emission between short GRBs 
and XRFs (X-ray Flashes, i.e. GRBs with E^;,so < 10^^ erg 
in Fig. Ill): in this picture, ?7^''°''' > 77;^^^. The two popu- 
lations are clearly distinct in terms of spectral peak energy 



during the prompt phase. 



with B^Jj""* 



> £, 



pk 



(Fig. ED. 



This may suggest that 77.^ anti-correlates with E-py^: this is 
further investigated in Sec. 14.1.21 

Short and long GRBs occupy different areas of the 
£x,iso vs. i?pk plane (Fig. El upper right panel) as well, 
demonstrating how the information from the X-ray LCs can 
be used to infer the GRB nature. Again, short GRBs fall 
below the long GRBs. 



4-0.2 The X-ray plateau and the prompt j-ray phase in 
long and short GRBs 

In the literature, the shallow decay (or "plateau") is as- 
sociated with an LC phase generally characterized by a 
mild slope (and a bsence of spectral evolution in the X-rays; 
iLiang et al.ll2007^ : this can be identified in type Ha and type 
HI light-curves. In type Ha (HI) GRBs, this phase starts 
at ti = ibi (ii = ib2) and ends at if = ib2 (if = ibs), 
with temporal slope a2 (as) and energy E2 (-E3). Short 
GRBs are under-represented in the class of GRBs showing 
clear evidence of plateaus in the X-rays. Only 2 short GRBs 
(out of 19 with C-like LCo 10%) possibly have plateaus: 
GRB051221A (Tgo = 1.4 s) and GRB070714B (T90 = 3 



^ 36 was the number of short GRBs in our starting sample; only 
19 of these have C-like LCs. 
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Table 4. From left to right: X and Y parameters to be correlated (the best fitting law reads: Log{Y) = q + mLog{X))\ best-fitting 
parameters as obtained accounting for the sample variance (D'Agostini 2005): slope (m), normalization (q), intrinsic scatter (tr); errors 
are given at 95% c.l. The last six columns list the value of the Spearman rank p, Kendall coefficient K and R-index r statistics and 
relative chance probability p. For each parameter couple, values reported in the first line refer to the entire sample, while in the second 
line we restrict our analysis to the long GRB class. X-ray fluences, fiuxes and observer frame times are computed in the 0.3-10 keV 
(observer frame) energy band; luminosities, energies and rest-frame times are computed in the 0.3-30 keV (rest-frame) energy band. 



X 


Y 


m 


9 


cr 


P 


P{p) 


K 


p(K) 


r 


p{r) 


Ey^Sso 


El^ 


1.07 ±0.03 


-4.1 ±72.2 


0.57 ±0.01 


0.70 


3 X 10-10 


0.52 


2 X 10-9 


0.79 


< 10-1" 






1.10 ±0.07 


-5.9 ±198 


0.58 ±0.01 


0.61 


3 X lO-'^ 


0.45 


6 X 10"'^ 


0.70 


10-9 


Sx.iso 


Lt 


1.21 ±0.06 


-15.6 ±169 


0.85 ±0.01 


0.58 


6 X lO-'^ 


0.45 


2 X 10"'^ 


0.70 


3 X 10-1" 






1.26 ±0.08 


-18.1 ±206 


0.85 ±0.01 


0.55 


3 X 10-" 


0.43 


8 X 10-'^ 


0.69 


10-9 


Sx,iso 


Li 


1.39 ±0.06 


-23.6 ±172 


0.83 ±0.01 


0.63 


3 X 10-8 


0.49 


2 X 10-8 


0.75 


< 10-1" 






1.37 ±0.08 


-22.8 ±212 


0.84 ± 0.02 


0.60 


3 X lO-'^ 


0.46 


2 X lO-'^ 


0.72 


10-1" 


tfF 


Lt 


-1.23 ±0.03 


51.9 ±0.46 


0.77 ±0.01 


-0.80 


< 10-"' 


-0.60 


< 10-1" 


-0.77 


< 10-1" 






-1.24 ±0.03 


52.0 ±0.45 


0.73 ±0.01 


-0.82 


< 10-1" 


-0.62 


< 10-1" 


-0.78 


< 10-1" 


Lt 


E2,X 


0.52 ±0.01 


26.8 ±11.8 


0.47 ±0.00 


0.67 


2 X 10-9 


0.51 


5 X 10-9 


0.80 


< 10-1" 






0.50 ±0.00 


27.6 ±10.8 


0.43 ± 0.00 


0.65 


10-8 


0.50 


2 X 10-8 


0.81 


< 10-1" 


E2,X 


El,K 


0.42 ± 0.02 


29.1 ±43.2 


0.81 ±0.01 


0.42 


4 X 10-"' 


0.29 


4 X 10-5 


0.45 


6 X 10-s 






— 


— 


— 


0.28 


6 X 10-3 


0.19 


0.06 


0.29 


6 X 10-3 


ti 


Ft 


-0.79 ±0.01 


-7.80 ±0.09 


0.45 ± 0.00 


-0.69 


< 10-1" 


-0.50 


< 10-1" 


-0.74 


< 10-1" 






-0.79 ±0.01 


-7.78 ±0.09 


0.45 ± 0.00 


-0.69 


< 10-"' 


-0.50 


< 10-1" 


-0.74 


< 10-1" 


-^7,iso 


Ex.iBO 


0.79 ±0.01 


10.0 ±20.6 


0.39 ±0.00 


0.86 


< 10-"' 


0.69 


< 10-1" 


0.88 


< 10-1" 






0.67 ±0.01 


16.5 ± 18.8 


0.29 ±0.00 


0.82 


< 10-1" 


0.63 


< 10-1" 


0.88 


< 10-1" 


-C'7,iso 


ttFL 
X 


0.89 ± 0.05 


3.85 ± 148 


0.65 ± 0.01 


0.64 


8 X 10-5 


0.48 


10-4 


0.74 


3 X 10-" 






0.93 ±0.10 


1.83 ± 287 


0.62 ± 0.02 


0.56 


10-3 


0.41 


10-3 


0.67 


6 X 10-5 


-^7,iso 


Si,x 


0.67 ±0.03 


15.9 ±91.3 


0.81 ± 0.02 


0.71 


2 X lO-'^ 


0.56 


2 X lO-'^ 


0.64 


5 X 10"" 






0.56 ±0.04 


21.6 ±126 


0.77 ±0.02 


0.62 


5 X lO--"' 


0.48 


4 X 10-5 


0.54 


5 X 10-4 


-C'7,iso 


E2,X 


0.92 ± 0.01 


2.96 ±33.5 


0.51 ±0.01 


0.76 


< 10-"' 


0.59 


< 10-1" 


0.85 


< 10-1" 






0.74 ±0.01 


12.6 ±35.1 


0.44 ± 0.00 


0.67 


10-8 


0.51 


2 X 10-8 


0.81 


< 10-1" 


-^7,iso 


Lt 


1.06 ±0.08 


-8.86 ±227 


1.03 ±0.04 


0.54 


9 X lO-"' 


0.41 


7 X 10-4 


0.70 


8 X 10-" 






1.05 ± 0.09 


-8.43 ± 260 


1.06 ±0.04 


0.50 


3 X 10-3 


0.37 


3 X 10-3 


0.68 


2 X 10-5 


^pk 


Ex,iBo 


0.98 ± 0.02 


49.5 ±0.15 


0.37 ±0.00 


0.63 


10-^ 


0.46 


5 X lO-'^ 


0.76 


< 10-1" 


-^pk 


Ex.iBO 


0.48 ±0.01 


27.0 ±16.4 


0.44 ± 0.00 


0.58 


2 X lO"'^ 


0.42 


3 X lO-'^ 


0.74 


< 10-1" 


^pk 


Lt 


0.86 ± 0.03 


2.35 ±87.6 


0.87 ±0.02 


0.50 


10-3 


0.39 


5 X 10-4 


0.76 


7x 10-8 


-^pk 


Ez.x 


0.60 ±0.01 


20.3 ±15.9 


0.43 ± 0.00 


0.58 


lo-'^ 


0.42 


3 X lO-'^ 


0.82 


< 10-1" 


S^ 


5x 


0.77 ±0.01 


-7.80 ±0.09 


0.45 ± 0.00 


0.79 


< 10-1" 


0.59 


< 10-1" 


0.77 


< 10-1" 






0.82 ± 0.00 


-1.58 ±0.10 


0.37 ±0.00 


0.78 


< 10-1" 


0.58 


< 10-1" 


0.78 


< 10-1" 


-^7,iso 


r llh 


0.71 ±0.01 


8.53 ±30.9 


0.55 ±0.01 


0.66 


3 X 10-1'' 


0.49 


1 X 10-9 


0.77 


< 10-1" 






0.54 ±0.01 


17.5 ±29.6 


0.45 ± 0.00 


0.55 


2 X 10-** 


0.40 


2 X 10-" 


0.70 


< 10-1" 


-L^7,iso 


J lOmin 
^X 


0.93 ±0.01 


-1.17 ±0.01 


0.45 ±0.01 


0.87 


< 10-1" 


0.67 


< 10-1" 


0.88 


< 10-1" 






0.78 ±0.01 


6.73 ±32.6 


0.40 ±0.00 


0.82 


< 10-1" 


0.63 


< 10-1" 


0.84 


< 10-1" 



s, extended emission not included). The corresponding per- 
centage for long GRBs is instead ~ 37%. 



4-0.3 The link between the X-ray luminosity and the 
prompt j-ray energy release 



The luminosity at the end of the plateau phase Lt is 
directly related to the total energy released in the second 
LC phase i?2,x (Table gj: E2,x oc Lf^^. It is interesting to 
note that of the two short GRBs, 070714B is a clear outlier, 
while 051221A is only barely consistent with the correlation. 
The peculiar GRB 060218 also shows a lower than expected 
£'2,x. iDainotti et a l. (2008) first reported a correlation be- 
tween Lt and tp^ for long GRBs. Here we confirm the cor- 
relation (with best-fitting Lt oc {tf'^)~^'^ Fig- and show 
that the two short GRBs with clear evidence of plateau are 
not consistent with the same scaling. 



The X-ray luminosity of the LC, Lx{t^^), correlates with 
the 7-ray energy released during the prompt emission for 
any t^^ between 100 s and 105 s. Here we arbitrarily se- 
lect two rest-frame times (Fig. |8} as an example. We find 
that the scatter of the correlation evolves with time, with 
the Lx{t^^) vs. Ej^iBo being tighter at early times (see Fig. 
[si. For this plot we require the GRBs to have been observed 
at those rest-frame times but relax the LC completeness 
requirement. No extrapolation of the observed LC is per- 
formed. 

At early times the LC luminosity tracks -E^.iso with 
limited dispersion around the best-fitting model Lx°™'" oc 
£'"fso. Short GRBs tend to lie below the best-fitting law of 
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Figure 6. Correlations involving E^^iso (1 — 10'' keV, rest-frame), -Bx,iso (0.3-30 keV, rest-frame) and the rest-frame prompt peak energy 
Ep]j. Black dots (red triangles): long (short) GRBs. Dashed line: best-fitting power-law model for the entire short plus long GRBs sample. 
Dotted line: best-fitting model for the long GRB class, only. The coloured and hatched areas mark the 68% confidence region around 
the best fit. Short GRBs and outliers are named. 
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Figure 8. 0.3-30 keV (rest-frame) X-ray LC luminosity measured at 10 minutes (left panel) and llhr (right panel) rest-frame vs. i?-,jso. 
Colour coding as in Fig. [5] The scatter of the correlation increases with time. 
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Log[tr/(s)] 

Figure 7. Luminosity at the end of the plateau phase, 0.3-30 keV 
(rest-frame) vs. end-time of the plateau. Colour coding as in Fig. 
[6] Blue points are for type la LCs. 



the long GRB class. When compared to the same relation 
at much later times (11 hours) we find that: (i) the rela- 
tion is now more scattered, suggesting that the X-ray LCs 



are more directly linked to the prompt 7-ray phase at early 
than at late times; (ii) while the relation is highly scattered, 
we note that all short GRBs of our sample lie below the 
long GRB relation: this is consistent with the steeper decay 
of the average short GRBs LC when compared to the long 
GRBs LC found in Sec 13.11 Our analysis therefore does no t 
confirm the previous results from iNvsewander et al.l l|2003 ). 
who found that short and long GRBs are consistent with the 
savne Lx vs. E^^iso scaling (note however that their Ej^iso is 
computed in a much narrower energy band). 



4-0.4 Observational biasses: temporal extrapolation 

The Swift re-pointing time ircp and end time of the observa- 
tions tend vary from GRB to GRB. Since Ex,iao is obtained 
by integrating the luminosity of each LC between trop and 
tend, one may wonder what is the effect of using different 
integration times for different GRBs. This is quantified as 
follows. To estimate the amount of energy lost at the end 
of the observations, we extrapolated the best-fitting profile 
of each GRB up to 10^ s (rest frame) and integrated the LC 
luminosity up until that time. Since GRBs m ay experience 
a jet break at late times (JRacusin et al.ll2009l ). this compu- 
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Figure 9. Impact of the temporal extrapolation of the observed 
LCs on the Ex iso vs. E^jgo correlation. Colour coding as in Fig. 
\E\ Light blue dots: i?x,iso values have been computed integrating 
the luminosity over a common rest-frame interval of time. 



tation may lead to an overestimate of the real energy lost. 
The amount of energy possibly lost at the beginning of the 
observationijj is estimated by conservatively extrapolating 
backwards in time the best-fitting profile to the minimum 
rest-frame Swift re-pointing time of our sample, which is 
12.5 s. For GRBs with T^^^ > 12.5 s, we adopt T^^ as the 
starting time for the integration to avoid extrapolating the 
luminosity to unrealistic values. This approach leads to an 
overestimate of the amount of energy lost before irep for 
the large majority of GRBs, as can be seen comparing the 
extrapolated temporal profile we are adopting here to the 
Swift-'BAT emission at the same rest-frame time (see e.g. the 
Swift Burst Anal yser BAT plus XRT LCs of GRB 050724, 
lEvans et al.ll20ld V The corrected -Ex.iso is shown in Fig. [5] 
(light blue points). Larger corrections (up to a factor ~ 9 for 
GRB 090510) are found to be applied to short GRBs. In spite 
of the very conservative approach we find that short GRBs 
are still either barely compatible or not consistent with the 
long GRB relation (as before), while the long GRBs relation 
is almost unaffected by this correction. 

We therefore conclude that in a logarithmic space the 
different rest-frame integration time used does not create 
or destroy correlations. The i?x,iso vs. -E^,iso correlation has 
been used here as an example: this result applies to all the 
relations presented in this paper. 



4.1 Multiparameter correlations 

We look here for correlations involving more than 2 param- 
eters (either from the X-rays or from the 7-rays). We first 
discuss the results from a Principal Component Analysis 
(PCA, Sect. l4.lTT|) and then show the existence of a tight 3- 
parameter correlation directly linking Exaso, iS^.iso and iJpk 
both in long and short GRBs (Sect. BT^ . 



* Note that the sample of C-like GRBs we use to look for corre- 
lations was pre-selected requiring an observed time of re-pointing 
tabs < 300 s to minimize this effect. 



Table 5. The three most significant PCs (85% of the total vari- 
ance) projected upon -E-^^isoi ^pki ^pki 'f's(f ^rid -Ex,iso- 





PCI 


PC2 


PC3 




40% 


66% 


85% 


^^,iso 


-0.561 


0.141 


0.171 


E^k 


-0.448 


-0.300 


-0.630 


Lpk 


-0.502 


-0.389 


- 


rpRF 
-'90 


-0.121 


0.794 


-0.512 


Ex,iso 


-0.466 


0.331 


0.588 



4-1.1 Principal Component Analysis 

The Principal Component Analysis (PCA) is a statistical 
technique designed to find patterns in data: it uses orthog- 
onal transformations to convert a set of possibly correlated 
variables into linearly uncorrelated (orthogonal) variables. 
Given a set of N events (GRBs in our case) described by M 
parameters, the PCA consists of the diagonalization of the 
covariance matrix: the eigenvenctors found are called prin- 
cipal components (PCs), while the eige nvalues consis t of the 
variance associated with each PC (see |jolliffell2002l . for de- 
tails). We performed a standardized PCA as recommended 
when the parameters have widely different variances: each 
parameter P is replaced by P = {Log[P] ~ Log[P])/aLog[p]- 
In this case the matrix to be diagonalized is not a covari- 
ance, but a correlation matrix. Calculations were performed 
using the statistical package RU. 

In Sect.Uwe showed that -Ex, iso is the X-ray parameter 
that still keeps information from the prompt 7-ray energy 
release. We now investigate its relation to other prompt pa- 
rameters, specifically -Epk, Lpk, Tgo^ (and -B^4so), using the 
PCA. This set of parameters is measured simultaneously in 
44 GRBs. Table [5] reports the three most significant PCs 
(86% of the total variance) projected upon the original 5 
variables. Each variable roughly contributes with compara- 
ble weight to the first PC; the second PC is instead domi- 
nated by fgo'". The third PC relates Exabo with Bpk. This 
result suggests that, while -Epk, ipk, Ex, iso and Ej^iso are in 
some way physically related to one another (see Sect. l4.lT2|) . 
the duration of the 7-ray energy release represents an addi- 
tional degree of freedom to the system. 



4-1.2 A GRB universal scaling: Ex, iso, i?7,iso and _Epk 

We look for a 3-parameter correlation involving -Ex, iso, E^^so 
and -Epk. The three variables are found to be correlated (see 
Fig. I10|) with the f ollowing best fitti ng law (obtained follow- 
ing the method bv lD'Agost"inil2005l ): 



Ex.ia 



10 



(0.58±0.25) 



-,(1.00±0.06) -, 



E, 



(0.60±0.10) 
pk 



(8) 



where -Ex i 



and Epk are in units of ergs and keV, 



respectively. The intrinsic scatter is as^ 
(1 a). We note that: 



0.30 ± 0.03 



This relation expands on the well known -Epk- -E^ 



http://www.r-project.org/ 
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Figure 10. Three parameter correlation involving -Bx,isoi ^7,iso a^nd -Ep^. Dashed line; best-fitting relation (Eq. [Sj; dot-dashed lines 
mark the 95% confidence area around the best-fitting law. Notably, long and short GRBs share the same scaling, with short and sub- 
energetic GRBs (like GRB 060218) occupying the same area of the plot. Inset: evolution of the efficiency parameter e = i?x,iso/^7,iso as 



a function of the spectral peak energy of the prompt emission Epi;. A reference e 
line. 



pk 



scaling has been marked with a blue dashed 



IjAmati et al.ll2008l and references therein) relation with the 
introduction of a third parameter (i?x,iso). 

• It combines information from the prompt and from the 
X-ray energy release which follows the prompt. While short 
GRBs are clear outliers of the -Epk- £'7,iso relation, they per- 
fectly fit into the -Bx,iso-i?7,iao--Bpk relation: the importance 
of the 3-parameter relation is that it combines short and 
long GRBs on a common scaling. As a result, considering 
the entire short plus long GRB sample, the scatter is re- 
duced by the introduction of the third variable (the intrinsic 
scatter of the Amati relation of our sample of long and short 
GRBs is (Je k = 0.37 to be compared to the intrinsic scatter 



of the 3-paramter relation on Spk which is cte 



: 0.29). 

Restricting our analysis to long GRBs, we find cte ^ = 0.17 
both for the -Epk- Ej^iao and for the 3-parameter correlation. 

• Short GRBs (like GRB051221A) and sub-energetic 
GRBs (like GRB 060218) occupy the same region of the 
-Epk- -E74so-Sx,iso space. The same is true for the peculiar 
long GRB 060614, la ter re-classified as a possible short GRB 
iGehrels et al.ll2006t ). In general, GRBs seem to divide into 
two groups with "normal" long GRBs occupying the upper- 
right area; short and peculiar GRBs together with XRFs 
(e.g. 050416A, 060218, 081007, 060614 also have a spectral 
peak energy below 60 keV) share the same lower-left region 
of the plot. 

• The best-fitting slope of the -Ex.iso vs. E^^iso relation 
of Fig. [6] reads: m = 0.79 ± 0.01 (see Table gj). The sig- 
nificant departure of m from 1 implies the more energetic 
long GRBs to have a lower e = -Bx,iso/i?7,iso oc l/E'^'J^oi 
with short GRBs being outliers of this relation. Interest- 
ingly, equation |8] implies: e oc y^tiu ; suggesting that the key 

parameter determining the 7-ray to X-ray ratio is not E~,^iso 
but the spectral peak energy _Bpk irrespective of the nature of 
the GRB (either long or short). This is clear from the inset 



of Fig. IIQI the higher the prompt peak energy the lower the 
e (th e GRB with E^k ~ 10"* keV is GRB 090510). We re- 
fer to IZhang et al.l ([20071) for a discussion of GRB radiative 
efficiencies derived from X-ray data. 

• The e(-Epk) scaling above can be interpreted as a phys- 
ical dependence of the radiative efficiency rj-y on Ep^: rj^ = 
Ej/{Ej -f Ek) ~ E-j/Ek oc E-t /Ej: as long as E ^ < Ek- 
This would imply 7/7 oc E°^ . See JFan et all i|2012l ) for a dis- 
cussion of this finding in the context of GRB photospheric 
models. Alternatively, a similar scaling could result for the 
long population if long GRBs with lower isotropic i?7,iso 
are less beamed than high energy GRBs during the prompt 
emission, but show otherwise similar beaming during the 
subsequent X-ray phase. In the first case, the e(i5pk) scal- 
ing would give direct information about the dissipative pro- 
cesses behind GRBs; in the second case, it would be an ob- 
servational effect, that nevertheless would provide valuable 
information about GRB jets and their opening angles. A 
complete and detailed discussion is beyond the scope of the 
present work and is provided by a companion paper (B12). 

• In Sec. 14.0.41 we showed that the different time intervals 
over which Ex,iso has been estimated do not severely affect 



the £'x,iso--E7 



E, 



pk. 



5 SUMMARY AND CONCLUSIONS 

We performed a comprehensive statistical analysis of Swift 
X-ray light-curves of 658 GRBs detected by XRT in the 
time period end of December 2004- end of December 2010. 
For the first time we present and analyse: (i) the properties 
of GRBs in a common rest-frame 0.3-30 keV energy band; 
(ii) we furthermore perform a comparative study of long 
and short GRBs; (iv) we cross-correlate the prompt 7-ray 
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properties and the X-ray LCs properties. We report below 
a summary of the major findings. 

From the spectral analysis of GRBs with redshift (Sec. 
[3:31): 



1. We find evidence for high intrinsic neutral hydrogen 
absorption NHhg ^ lO^^cm"^ even at 2: < 2. The average 



10^ 



cm (mean value of 



value for long GRBs is NHhg 
the logarithm of the NHhg). 

2. Short GRBs map the low end of the distribution with 



mean NHhg 



10^ 



cm . However, there is no evidence 



for short GRBs to show a lower NHhg when compared to 
long GRBs in the same redshift bin. 

The analysis of 297 long GRBs with complete X-ray 
light-curve^^ reveals that: 

3. The average energy released in X-rays (0.3-30 keV, 
rest-frame) is Ex,iBo ~ 7 x 10^^ erg typically representing 
~ 7% of Ej^iso (Sec. 13. 2p . The two quantities are statisti- 
cally correlated: -Ex iso oc E^'f^^. Also: -Ex iso oc -Epk (Sec. 

mm . 

4. The -Ex, iso distribution does not extend beyond 10^^ 
erg (Sec. 13. 2[) possibly suggesting the existence of a 
maximum available energy budget (the record holder is 
GRB 080721 with Ex.iso ~ 10^^ erg). Also, for z > 2 we 
are not sensitive to the population of GRBs with -Ex, iso < 
lO^'^ erg, so that the low-energy tail of the -Ex, iso distribution 
is currently under-sampled. 

5. The X-ray luminosity of the LCs at any rest frame 
time between 100 s and 10^ s is found to correlate with 
E-y^iao (Sec I4.0.3|) : the scatter of this correlation increases 
with time which might suggest that early time X-rays are 
more tightly related to the prompt phase. 

In the case of short GRBs (19 have C-like LCs): 

6. The median luminosity light-curve of short GRBs (Sec. 
13. ip is a factor ~ 10 — 30 dimmer than long GRBs in the rest- 
frame time interval 10^ — 10* s, has a steeper average decay 
(oc t~^'^ vs. oc t~^) and shows no evidence for clustering at 
late times (contrary to long GRBs). 

7. Short GRBs populate the low-energy tail of the Ex, iso 
distribution, with E^°^^ ~ ■?7)^x?so ^^'^ ^^ average 
£^x;'i°so' ~ 10^° erg (Sec. 1321). Short GRBs are more energy 
deficient during the second LC phase when compared to long 
GRBs. 

8. Short bursts are clear outliers of the Ex, iso — E-,^iso and 
£x,iso ~ -Epk relations established by the long population, 
with E^°^J a factor > 50 below expectations (Sec. I4.0.ip . 
Short GRBs are also found to lie below the L^^^ vs. -E^,iso 
relation established by the long class. 

9. Short GRBs are under-represented in the class of GRBs 
showing clear evidence of plateaus in the X-rays. Only 2 
GRBs out of 19 possibly have plateaus (10%). The corre- 
sponding percentage for long GRBs is instead ~ 37%. While 
the limited sample size does not allow us to draw firm con- 
clusions, we note that X-ray plateaus are more commonly 
detected in long GRBs (Sec. 14.0. 2p . 

10. The two short GRBs with X-ray LC plateaus in our 
sample are outliers of the Lf vs. tf^ relation fSec. l4.0T2|) . 



Irrespective of the long or short GRB nature, we find no 
statistically significant correlation involving the rest frame 
prompt duration Tc^^, the intrinsic column density NHhg 
or the temporal slopes of the X-ray LCs (Sec.[4|. The Tg^^ 
basically accounts for the second strongest Principal Com- 
ponent (Sec. l4.lT]) . suggesting that while -Epk, -Lpk, E^x,iso 
and E^^iso are related to one another, the 7-ray duration 
represents an additional degree of freedom to the system. 

We showed in Sec. 14.1.21 the existence of a 3-parameter 
correlation that links Ex, iso, -E^.iso and -Epk: -Ex, iso oc 

^pk 



/pi. 00 /p0.60^ 

\ 7,iso/ -^p 



(i) Short and long GRBs share the same scaling. 



„^''°° oc E_k which can be 



(ii) This correlation implies 

interpreted as rjj oc -Epjf (where rij is the radiative effi- 
ciency) . 

(iii) Standard long GRBs and short GRBs (together with 
peculiar GRBs and XRFs) occupy a different region of the 
E?x,iso-E?7,iso--Epk plane. 

The results from our analysis are publicly available r'^ I 
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APPENDIX A: GLOSSARY 

This section provides the list of symbols used. As a gen- 
eral note: X-ray energies (fluences) were computed from the 
time of the Swift-XI^T repointing up until the end of the 
observation; no temporal extrapolation was performed. The 
values reported assume isotropic emission. X-ray fluences 
and fluxes are reported in the 0.3-10 keV (observer-frame) 
energy band; energies, luminosities and intrinsic time scales 
are computed in the 0.3-30 keV (rest-frame) band. 



• an'- temporal slope of the normal decay phase. Type la: 
Qn = 02; type Ila: On = as; type III: On = 04. 

• cvst: temporal slope of the steep decay phase. Type lb 
and Ila : a^t ~ cti; type lib and III: ast ~ oti (see Fig. [l|. 
The zero-time of the power-law decay is assumed to be the 
BAT trigger time (i. e. to = 0). 

• oi^x''^ : temporal slope of the steep decay phase assuming 
io ~ Tgo- 

• ftsh: temporal slope of the shallow decay (or plateau) 
phase. This corresponds to 02 and 0^3 for type Ila and type 
III light-curves, respectively. 

• Fa,: XRT 0.3-10 keV (observer frame) spectral photon 
index from this paper. 

• E-, 



prompt emission in the rest-frame 1 — 10 keV energy band 



v7,iso. isotropic equivalent energy released during the 
pt emission in the r e 
from lAmati et al.l (|2008f l. 

• i?pk: rest-frame peak energy of the vFy spectrum during 
the prompt 7-ray emission from lAmati et al.l l|2008h . 

• Ff (Lf): flux (luminosity) at the end of the plateau (i.e. 
at t = tf). 

• Fi (Li): flux (luminosity) at the beginning of the plateau 
(i.e. at i = ii). 

• ipk.iso: 1 — 10^ keV (rest frame) i sotropic peak lum i- 
nosity during the prompt emission from lNava et al.l (|2008l ) . 

• L^x^: luminosity at 11 hours rest-frame. 

• 1/3?°""- luminosity at 10 min rest-frame. 

• NHtot: total neutral hydrogen column density. 

• NHhg: intrinsic neutral hydrogen column density at the 
redshift of the GRB. 

• S'i,x (-Ei,x): fluence (energy) released during the first 
phase of the X-ray light-curve. Type lb and Ila: i?i,x = E\\ 
type lib and III: -Ei,x = -Ei 4- -E2. Fluences follow the same 
definition scheme. E\, E2, -E3 and E4 has been defined fol- 
lowing Fig. [1] 

• 52, x (F2,x): fiuence (energy) released during the second 
phase of the X-ray light-curve. Type la: -E2,x ~ Ei + E2; 
type lb: E2,x = E2; type Ila: _E2,x = E2 + E3; type lib: 
-E2,x = E3; type III: £;2,x = E2 + E4 (see Fig. [TJ. Same 
definition scheme for fiuences. 

• S-y (£^^"^^°): 15-150 keV (observer frame) fluence (en- 
ergy) released during the prompt emission as calculated by 
ISakamoto et all l|201ll ). 

• 5x (Fx.iao): X-ray fluence (energy). 

• Sx^ (E^): X-ray fluence (energy) associated to flares. 
For each GRB, the total fluence (energy) released in X-rays 
reads: S^'^+Sx (F^^^+Fx.iso). 

• Sf^, SEk {Efk, ^2,x): X-ray fluence (energy) of flares 
superimposed on the first and second light-curve phase. 

• if, t^^ , t™°: end time of the plateau phase: observer 
frame, rest frame and in T90 units. This parameter corre- 
sponds to i(,2 and tbs for type Ila and type III light-curve, 
respectively. 



ti, if 



ii 



start time of the plateau phase: observer 



frame, rest frame and in T90 units. This parameter corre- 
sponds to i(,i and t(,2 for type Ila and type III light-curve, 
respectively. 

• At: plateau duration defined as tf — ti. 

• Tgo , T^ : duration of t he 15-150 keV prompt emission 
from lSakamoto et al.l l|2011f) . in the observer and in the rest- 
frame, respectively. 
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Table Bl. Best-fitting parameters of the 0.3-10 keV (observer frame) light-curves in de-absorbed flux units, as obtained following the 
procedure outlined in Sec. 12. ll GRB name, LC type (as defined in Sec. 12. it . redshift, power-law indices (oi, 02, 03, 04) and errors, 



n i201ll : 



to lSakamoto et al.ll201ll for GRBs detected before December 2009, and to the refined BAT GCNs otherwise), power-law index of the first 
segment when to = '^90 (0^1^^") a^nd error, x^, degrees of freedom, p-value. Normalizations are given in 10"^'' erg cm"-^ s~^, times in 
seconds. A redshift equal to indicates that no reliable estimate of this parameter is available from the literature. For the other columns, 
—9 indicates that the value is absent (i.e. there is no such LC phase). Note on the LC type column: C-GRBs (i.e. GRBs with complete 
LCs) are defined as promptly repointed GRBs t^ep < 300s whose fading was followed up to a factor 5-10 from the background limit; if 
this is not the case, the GRB is flagged as U-like (i.e. GRB with truncated LC). The flag F (N) indicates that flares have (have not) been 
detected. This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance. 



GRB 


Type 


z 


ai 


a"ai 


a2 


a"a2 


Q3 


0"a3 


04 


rjdi 


Log[ti,i] 


ctfci 


Log\ti,i\ 


ft 62 


iog[t(,3] 


"■(bS 


041223 


OUN 





1.91 


0.44 


-9. 


-9. 


-9. 


-9. 


-9 


-9 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


050124 


OUN 





1.44 


0.17 


-9. 


-9. 


-9. 


-9. 


-9 


-9 


-9. 


-9. 


-9. 


-9. 


-9 




-9. 


050126 


ICN 


1.29 


2.5 


0.51 


0.862 


0.25 


-9. 


-9. 


-9 


-9 


2.72 


0.4 


-9. 


-9. 


-9 




-9. 


050128 


lUN 





0.758 


0.12 


1.38 


0.14 


-9. 


-9. 


-9 


-9 


3.69 


0.38 


-9. 


-9. 


-9 




-9. 


050219A 


ICN 





3.68 


0.36 


0.779 


0.11 


-9. 


-9. 


-9 


-9 


2.4 


0.054 


-9. 


-9. 


-9 




-9. 


050219B 


OUN 





1.42 


0.044 


-9. 


-9. 


-9. 


-9. 


-9 


-9 


-9. 


-9. 


-9. 


-9. 


-9 




-9. 


050315 


IICN 


1.95 


-0.295 


0.18 


0.819 


0.037 


3.83 


0.18 


-9 


-9 


2.71 





3.82 


0.14 


-9 




-9. 



Table B2. Continued from Table iBll 



Nx 


"■jvi 


N-2 


0'JV2 


SI 


S2 


Tbo 


„T90 


O"aT90 


2 

X 


d.o.f. 


p-val 


1.45e7 


6.3e7 


-9. 


-9. 


-9. 


-9. 


109. 


1.91 


0.44 


8.39 


25. 


0.999 


4.72e4 


8.e4 


-9. 


-9. 


-9. 


-9. 


3.93 


1.44 


0.17 


31.6 


30. 


0.384 


0.0673 


0.11 


-9. 


-9. 


0.5 


-9. 


48. 


2.5 


0.51 


3.54 


10. 


0.966 


0.458 


0.49 


-9. 


-9. 


-0.5 


-9. 


28. 


0.758 


0.12 


69.8 


148. 


1. 


0.457 


0.14 


-9. 


-9. 


1. 


-9. 


23.8 


3.13 


0.31 


56.6 


79. 


0.973 


1.26e5 


4.8e4 


-9. 


-9. 


-9. 


-9. 


28.7 


1.42 


0.044 


81.9 


127. 


0.999 


0.179 


0.027 


2.15e9 


1.8e9 


-0.5 


-9. 


95.6 


-0.295 


0.18 


155. 


218. 


1. 



Table B3. 0.3-10 keV (observer-frame) fluence table. From left to right: GRB name, LC type (as defined in Sec. 12. it . redshift, initial 
^min) and final {Tmax) time of the observations, total fiuence {Sx,iso) with error, fluence of the different LC phases (Si, S2, S3, S4) 
and errors, fluence of the flares in different parts of the LC (S^^, 82^, S^^, S^^) and errors. Fluences are given in erg cm~^. A 
redshift equal to indicates that no reliable estimate of this parameter is available from the literature. A "—9" indicates that the LC 
does not show such phase and the value of that parameter is therefore absent. Finally, for the columns containing information from flares 
superimposed on the power-law decay, indicates that no statistically significant positive fluctuation has been detected. This table is 
available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance. 



GRB 


Type 


z 


Log[T.,ni„] 


Log[T„,ax] 


Sx,^so 


<^sx 


Si 


""SI 


52 


0-S2 


53 


<^S3 


S4 


0-S4 


041223 


OUN 





4.22 


4.45 


8.8e-8 


6.7e-9 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


050124 


OUN 





4.05 


6.4 


1.6e - 7 


3.4e - 8 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


050126 


ICN 


1.29 


2.12 


4.83 


3.7e-8 


6.8e - 9 


1.5e-8 


2.e-9 


2.2e-8 


6.3e - 9 


-9. 


-9. 


-9. 


-9. 


050128 


lUN 





2.23 


4.84 


7.9e - 7 


2.4e-8 


4.6e - 7 


1.8e-8 


3.3e-7 


1.5e-8 


-9. 


-9. 


-9. 


-9. 


050219A 


lUN 





2.05 


6.2 


3.5e-7 


1.4e-7 


l.le -8 


2.e - 9 


3.1e-7 


l.le - 7 


-9. 


-9. 


-9. 


-9. 


050219S 


OUN 





3.5 


6.22 


9.4e - 7 


6.1e-8 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


050315 


IICN 


1.95 


1.92 


5.93 


l.le -6 


3.3e - 8 


2.4e - 7 


1.4e - 8 


7.e-8 


8.8e - 9 


8.e-7 


2.8e-8 


-9. 


-9. 



Table B4. Continued from Table lB3l 



Sf^ 



TgiFL 



Sf^ 



^S2P^ 



si^ 



sf^ 



crs4.Fi. 



-9. 
-9. 






-9. 





-9. 
-9. 






-9. 





-9. 


-9. 


-9 


-9 


-9. 


-9. 


-9 


-9 


-9. 


-9. 


-9 


-9 


-9. 


-9. 


-9 


-9 


-9. 


-9. 


-9 


-9 


-9. 


-9. 


-9 


-9 








-9 


-9 
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Table B5. 0.3-30 keV (rest-frame) energy table. From left to right: GRB name, LC type (as defined in Sec. l2.l(l . redsliift, initial (Tmin) 
and final (Tmax) time of the observations, total energy {Ex iso) with error, energy of the different LC phases (-Bi, E2, -E3, -B4) and 
errors, energy of the flares in different parts of the LC (E^^, -Bj^^, E^^ , E^^) and errors. Energies are given in erg. A redshift equal 
to indicates that no reliable estimate of this parameter is available from the literature. A "—9" indicates that the LC does not show 
such phase and the value of that parameter is therefore absent. Finally, for the columns containing information from flares superimposed 
on the power-law decay, indicates that no statistically significant positive fluctuation has been detected. This table is available in its 
entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content. 



GRB 


Type 


z 


LoglTm.n] 


iog[T,„„x.] 


£x,..o 


ffiX 


E, 


0"B1 


£2 


""£2 


E3 


<^E3 


Ei 


"■£4 


050126 


ICN 


1.29 


2.12 


4.83 


3.7e - 8 


6.8e-9 


1.5e-8 


2.e-9 


2.2e-8 


6.3e - 9 


-9. 


-9. 


-9. 


-9. 


050315 


IICN 


1.95 


1.92 


5.93 


l.le - 6 


3.3e-8 


2.4e - 7 


1.4e-8 


7.6-8 


8.86 - 9 


8.6-7 


2.8e-8 


-9. 


-9. 


050318 


OUN 


1.44 


3.52 


5.65 


2.5e-7 


l.le -8 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


050319 


IICN 


3.24 


1.98 


6.15 


5.e-7 


5.6-8 


2.5e-8 


3.6e-9 


1.7e-7 


6.76 - 9 


3.6-7 


4.9e - 8 


-9. 


-9. 


050401 


ICF 


2.9 


2.14 


5.9 


1.2e-6 


6.2e-8 


5.e-7 


7.5e-9 


7.3e - 7 


6.16-8 


-9. 


-9. 


-9. 


-9. 


050408 


OUN 


1.24 


3.41 


6.47 


5.4e - 7 


3.9e-8 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 


-9. 



Table B6. Continued from Table iBSl 



Bf^ 


"eifl 


Bf^ 


CTE2Fi 


E^^ 


0-£;3FL 


£f^ 


0-E4F1, 














-9. 


-9. 


-9. 


-9. 




















-9. 


-9. 








-9. 


-9. 


-9. 


-9. 


-9. 


-9. 




















-9. 


-9. 














-9. 


-9. 


-9. 


-9. 








-9. 


-9. 


-9. 


-9. 


-9. 


-9. 
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